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Linear  Envelopes  for 
Uniform  B-spline  Curves 

David  Lutterkort  and  Jorg  Peters 


Abstract.  We  derive  an  efficiently  computable,  tight  bound  on  the 
distance  between  a uniform  spline  and  its  B-Spline  control  polygon  in 
terms  of  the  second  differences  of  the  control  points.  The  bound  yields 
a piecewise  linear  envelope  enclosing  the  spline  and  its  control  polygon. 
For  quadratic  and  cubic  splines  the  envelope  has  minimal  possible  width 
at  the  break  points,  and  for  all  degrees  the  maximal  width  shrinks  by  a 
factor  of  4 under  uniform  refinement.  We  extend  the  construction  to  tight 
envelopes  for  parametric  curves. 


§1.  Motivation  and  Overview 

The  central  feature  that  allows  reasoning  about  nonlinear  piecewise  polynomi- 
als is  the  fact  that  a spline  is  closely  outlined  by  its  B-spline  control  polygon. 
The  efficiency  of  many  applications  depends  crucially  on  a good  estimate  of 
the  distance  separating  spline  and  control  polygon.  For  rendering,  a refined 
control  polygon  is  rendered  instead  of  the  curve  itself.  For  curve-intersection 
an  efficient  and  robust  technique  is  to  recursively  refine  and  intersect  control 
polygons  [2].  Assessing  the  exactness  of  these  operations  requires  a uniform, 
linear  bound  on  the  distance  of  the  curve  and  its  (refined  B-spline)  control 
polygon.  The  efficiency  is  improved  if  the  effect  of  the  refinement  can  be  pre- 
dicted rather  than  just  measured.  Of  the  two  classical  bounding  contructs, 
axis-aligned  min-max  coefficient  boxes  and  the  convex  hull,  the  first  yields 
only  a loose  envelope  and  neither  yields  a priori  estimates. 

This  paper  introduces  quantitative  bounds  that  can  be  computed  more 
efficiently  than  convex  hulls,  and  yield  a simple  piecewise  linear  envelope  en- 
closing spline  and  control  polygon  (see  Figure  1)  whose  maximal  width  con- 
tracts to  1 / 4th  when  the  knot  spacing  is  halved.  The  computation  of  the 
envelope  of  a degree  d spline  consists  of  computing  the  second  differences  of 
its  control  points  and  looking  up  or  calculating  d — 1 constants,  the  values 
of  a fixed  set  of  splines.  The  sum  of  the  constants,  ( d + l)/24,  provides  a 
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Fig.  1.  A cubic  curve  (black)  and  its  control  points  (black  squares).  The  enve- 
lope (grey)  is  constructed  with  the  bound  from  Theorem  2. 

second,  even  simpler,  but  generally  much  coarser  bound  (Figure  2).  Both 
bounds  are  piecewise  linear  with  breaks  at  the  corners  of  the  control  polygon 
and  are  sharp  for  quadratic  and  cubic  splines  in  the  sense  that  at  every  corner 
of  the  control  polygon  the  distance  between  the  spline  and  polygon  is  matched 
exactly. 

This  paper  derives  these  bounds  for  functions  and  establishes  the  conver- 
gence of  the  bound  under  uniform  refinement;  the  bounds  are  then  applied  to 
curves  to  obtain  localized  envelopes. 

§2.  Notation 

A scalar-valued  piecewise  polynomial  p of  degree  d is  a uniform  B-spline  if  it 
can  be  represented  as 

p=J2bkNk,  6*  el R,  Nk  = N(--k), 

fee  2z 

where  N is  the  B-spline  of  degree  d supported  on  the  interval  [0,  d + 1)  and 
with  the  uniform  knot  sequence  TL  (c.f.  [2]).  For  simplicity,  we  assume  that 
both  the  control  point  sequence  and  the  knot  sequence  are  biinfinite. 

The  control  polygon  t of  p is  the  piecewise  linear  interpolant  of  the  control 
points  bk  at  the  Greville  abscissae 

tk  — k + (d  + l)/2. 

Over  the  interval  the  control  polygon  is  (.(t)  = Lk(t ; bk ,bk+l)  where 

we  denote  the  line  segment  from  (t£,ai)  to  (f£+1 , 02)  by 

Lfc(t  ;ax,a2)  = ai(<k+i  - t)  + a2{t  - t*k). 

The  linear  interpolant  of  a function  / over  this  interval  will  be  abbreviated 
as  Lfc  (/)  = bk(- f{tl),  f(tl+1)).  The  (centered)  second  differences  of  b are 
defined  as 

A26’  = 6i_1  - 26*  + 6i+1. 

The  first  and  the  last  basis  function  that  are  supported  on  K.tfc+il  are  N~ 
and  Nk  with 

k = k + l-[d/2\,  k = k-l+[d/2\. 
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Fig.  2.  A cubic  curve  (black)  and  its  control  points  (black  squares).  On  top 
the  envelope  (grey)  is  constructed  with  the  bound  from  Theorem  3,  on 
the  bottom  from  the  tighter  bounds  of  Theorem  2. 


§3.  Uniform  B— splines  Bounds 

The  key  observation  for  deriving  the  bounds  is  that  the  difference  between  a 
uniform  B-spline  p and  its  control  polygon  can  be  factored  into  two  parts:  the 
second  differences  of  the  control  polygon  and  splines  pki , which  are  indepen- 
dent of  p. 

Theorem  1.  Over  the  interval  [t£,tfc+1],  the  difference  between  a uniform 
B-spline  p and  its  control  polygon  l is  given  by 


p-e  = ]TA2bi/3ki, 


Pki  = 


\Z7=i(j-i)Nj  i>  k. 


The  functions  Pki  are  non-negative  and  convex  on  the  interval  [f£,  t£+1]  and 
Pki(t*k ) > 0 if  and  only  if  i € [k,  k]. 

Proof:  We  write  p — i over  [t£, <£+1J  as 


k k 

E bi  = E W)  - M* ; Sik,Si,k+i)) , 

i=k  i=k 

where  6ik  = 1 if  i = k and  0 otherwise.  We  show  that  aki  = A2Pki  '■  the 
partition  of  unity  ]TT  Nl  = 1 implies  that  aki  = 0 and  the  linear  precision 
of  B Splines,  Xw  t* -W’(i)  = t implies  on  the  interval  [t*kP*k+i\  that  J2i  i-aki  — 0. 
Hence,  for  any  i,  ~ i)akj  = 0. 

For  i > k, 

k k i 

Pki  = E(j  - i)Nj  - E(i  - *)«*i  = E(*  ~ 

j=i  j=k  j=k 
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so  that  Pki  = — j)akj  for  any  i.  It  is  now  straightforward  to  verify 

that  A2 Pki  — &ki  and  summation  by  parts  yields 

k k k 

p-ek  = YJbiaki  = YJbi  ^Pki  = E /?«■ 

i=Jk  *=fc  »=fc 

The  functions  (3ki  are  non-negative  since  their  B-spline  coefficients  are  non- 
negative. The  convexity  of  the  f3ki  over  [t*k,  t£+1]  follows  from  the  convexity  of 
their  B-spline  control  polygons:  for  i > k,  the  part  of  the  control  polygon  of 
[3ki  that  influences  (iki  over  [fjjl,  ] lies  on  the  function  max{-  — f*,0}  while 
for  i < k it  lies  on  max{£*  — •,  0}.  In  both  cases,  the  control  polygon  of  /?*,,, 
and  hence  pki,  is  non-negative  and  convex.  □ 

Theorem  1 immediately  gives  us  a piecewise  linear  envelope  on  p — t. 

Theorem  2.  Over  the  interval  [fj,  the  difference  between  a uniform 

B-spline  p and  its  control  polygon  f is  bounded  by 

Lk  (]T  A 2v  &,)  <p-e<Lk(22  A }b*  Phi)  , 

where  A 2b'  = max{A2fc\0}  and  A2b'  = min{A2b',0}. 

Proof:  We  have  from  Theorem  1 

P - e = E A^  fai  = E A2  + E A2  6< 

i i t 

The  positivity  of  the  pki  implies  that  the  first  sum  on  the  right-hand  side  is 
positive  and  the  second  is  negative  and  therefore 

YJ^biPki<p-^<Y,AtbiPki- 

i i 

Since  the  /3ki  are  convex  over  they  can  be  bounded  linearly  to  yield 

the  bound  of  Theorem  2.  □ 

An  even  simpler  envelope  can  be  derived  by  bounding  the  sum  of  the  pki 
at  t*k  by  the  constant  ( d + l)/24. 

Theorem  3.  Over  the  interval  [t£, the  difference  between  a uniform 
B-spline  p and  its  control  polygon  i is  bounded  by 

\p-e\<^Lk(-;\lA2bHk,\\A2b\\k+1), 

where  ||  A26H jt  = max{|A2&'|  : i € [k,  fc]}.  If  d = 2 or  d = 3 the  relation  holds 
with  equality  at  the  t*k. 

Proof:  By  Theorem  1 and  the  convexity  of  the  (3ki  over  [t£,  f£+1],  we  have 
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Fig.  3.  A self-intersecting  quartic  curve  (black)  and  its  control  points  (black 
squares).  The  envelopes  (grey)  are  constructed  with  the  bound  from 
Theorem  2.  The  envelope  converges  rapidly  to  the  curve  as  the  com- 
parison of  the  original  envelope,  top,  and  the  envelope  after  one  step  of 
uniform  refinement,  bottom,  shows. 

The  theorem  follows  if  we  can  show  that  JT  Ski  < [d  + l)/24: 

1 Z&ki  = Y ]0*  ~ j)  ak>  = Y Y^  “ j)  a*b  = Y aki  Y * 

i i=k  j—k  j=k  i=j  j i— 0 

j ' 7 jw 

Regardless  of  the  degree  of  p,  z is  the  quadratic  polynomial 


Since  z is  a positive  and  convex  function,  z attains  its  maximum  over 
at  one  of  the  endpoints  of  the  interval.  Its  values  there  are 

*<&)  = *&+ i)  = ^r 

and  hence  z(t)  = Ski  < (d  + 1) /24  for  all  t G [fjt,  ^+il- 

The  number  of  Ski  that  are  nonzero  at  t\  is  d — 1 for  d even  and  d — 2 
for  d odd,  i.e.  only  Skk  is  nonzero  at  t*k  if  p is  quadratic  or  cubic.  But  then  all 
inequalities  of  equation  (1)  become  equalities  as  claimed.  □ 

Computing  the  bounds 

To  compute  the  bounds  for  quadratics  or  cubics  no  B-spline  evaluation  is 
required,  since  only  Skk(t*k)  = (d  + 1) /24  is  nonzero.  For  d > 3,  it  suffices  to 
look  up  tabulated  values  Ski{t*k)  f°r  ~\A/2\  < i < [d/2\ . Forming  the  inner 
products  of  Theorem  2 and  Theorem  3 at  t*k  is  straightforward. 
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§4.  Uniform  Refinement 

An  important  operation  on  B-splines  is  the  refinement  of  the  knot  sequence 
or  knot  insertion.  Knot  insertion  changes  the  representation  of  the  piecewise 
polynomial  p over  the  original  knot  sequence  to  one  over  a larger  knot  sequence 
and  reduces  the  distance  between  spline  and  control  polygon  (c.f.  Figure  3). 

After  halving  the  distance  between  knots  the  new  control  points  bk  in 
p(t)  = J2k  bkNk(t)  = E/t  bkNk(2t)  are  given  by 


Theorem  4.  The  second  differences  A26'  of  the  refined  control  polygon  are 
bounded  by  the  second  differences  A2 b‘  of  the  original  control  polygon 

max  |A26'|  = - max  |A26'|. 
i 4 » 

Proof:  The  second  derivative  p"  of  p is  given  by 

p"(t)  = J2^bk~1Nk_2(t)  = £A  2bk~lNkd_2k{2t), 

k k 


which  means  that  the  A26*  can  be  obtained  from  the  A26‘  via  (2)  as 
2dA262'  = £ A26‘->,  2dA262<+1  = £ (d“  ^A^"'. 


The  proof  follows  from  Ej  (*/)  = 2d  1 and  T,j  ('2/)  = Ej  (2*_\)  = 2<i  2 • 

□ 

Theorem  4 yields  the  following  a priori  estimate  on  the  number  of  subdi- 
visions er  needed  to  bring  spline  and  control  polygon  within  a given  distance  e: 


o(p,s)  = [log4 


(d+l)||A2frl| 
24s  1 


Examples:  For  quadratic  B-Splines,  uniform  refinement  is  called  Chaikin’s 
algorithm,  and 

b2i  = 2_2(36'_1  + V),  b2i+1  = 2_2(6’_1  + 36;). 


This  yields 

A262i  = A262i_1  = ^A2&i_1, 

4 

i.e.  every  second  difference  is  guaranteed  to  decrease  by  a factor  of  four.  Sim- 
ilarly, for  cubic  B-Splines  we  have 

b2i  = 2~3(6‘~2  + 66' _1  + 6'),  62i+1  = 2_3(46'_1  + 46*), 


and 


A262i  = -A26'_1, 


A262i+1 


1 + A26’ 

4 


2 
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Fig.  4.  Constructing  the  envelope  of  a curve  from  the  bounding  rectangles  Sk 
and  Sk+1\  only  the  outer  line  segments  tij  and  u\  are  part  of  the  convex 
hull  of  Sk  U Sk+1  and  the  envelope. 

§5.  Curve  Envelopes 

A parametric  curve  p is  in  uniform  B-spline  form  if  p = ]T\-  b3  N J where  the 
fr7  € 1R”  are  the  control  points  of  p and  the  uniform  B-spline  basis  N3  is 
defined  as  in  Section  2.  The  curve  p is  closed  if  the  control  point  sequence 
( b ■’)  is  periodic. 

The  functional  bounds  are  applied  componentwise  to  parametric  curves. 
Then  each  control  point  and  the  curve  point  corresponding  to  its  Greville 
abscissa  lie  in  a box  whose  width  in  the  ith  component  is  the  bound  in  the 
ith  component.  It  is  now  convenient  to  restate  the  bounds  from  Theorems  2 
and  3 more  abstractly  as 

e(t)  < p{t)  - £{t)  < e(t)  for  t € [t*k,t*k+ 1].  (3) 

For  curves  p,  the  bound  in  the  i-th  component  is  denoted  by  e{  <Pi~li  < e*. 
By  (3),  p{tk)  is  located  in  the  axis-aligned  box  Sk, 

Sk  = {x  | ej(tfe)  <Xi-bk  < ei(t*k ) for  all*  = 1, . . . ,n}. 

Each  point  of  the  curve  segment  p(t),  t € \tk,t*k+ J,  lies  in  a box  S{t), 
that  by  the  linearity  of  e and  e is  a convex  combination  of  Sk  and  Sk+1: 

S(t)=Lk(t-,Sk,Sk+1). 

The  curve  segment  is  therefore  contained  in  the  union  of  all  S(t),  t € [t*k, 
which  is  the  convex  hull  Hk  of  the  corners  of  Sk  and  Sk+1.  To  be  specific,  we 
discuss  the  case  of  planar  curves. 

Enveloping  planar  curves 

Let  vk,i=  1, ...  ,4,  be  the  line  segments  connecting  corresponding  corners  of 
Sk  and  Sk+l;  that  means  vk  connects  the  lower  left  corner  of  Sk  to  the  lower 
left  corner  of  Sk+1,  vk  connects  the  lower  right  corner  of  Sk  to  the  lower  right 
corner  of  Sk+1  etc.  as  in  Figure  4. 

Hk  consists  of  parts  of  the  boundaries  of  Sk  and  Sk+1  and  exactly  two 
additional  line  segments  uk  and  it 2 chosen  from  the  vk . Since  uk  and  112  are 
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part  of  the  convex  hull  Hk,  they  do  not  intersect  the  interiors  of  Sk  and  Sfc+1 . 
We  do  not  need  to  actually  compute  intersections  of  the  vk  and  Sk,  Sk+l  to 
select  uk  and  uk:  since  Sk  and  are  axis-aligned  it  suffices  to  look  at  the 
signs  of  the  slopes  of  the  vk . The  uk  are  separated  by  the  line  from  bk  to 
bk+1 ; we  call  the  one  lying  to  the  left  of  this  line  uk  and  the  one  lying  to  the 
right  of  this  line  uk. 

The  sets  Ui  — {u^}  are  not  yet  polylines:  consecutive  line  segments  uk 
and  uk+1  may  intersect  or  not  touch  at  all.  But  note  that  the  line  extending 
uk  always  intersects  the  one  extending  u*-1.  We  obtain  a proper  polyline 
Wi  with  exactly  one  line  segment  for  each  control  point  of  p by  taking  this 
intersection  as  starting  point  and  the  intersection  with  the  line  through  uk+1 
as  the  end  point  of  W{.  The  polylines  Wi  and  W2  then  form  a local  envelope 
of  p:  the  curve-piece  p([t^,t^+i])  lies  in  the  quadrangle  spanned  by  the  k- th 
pieces  of  Wi  and  W2 . 
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